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When the dynamics of a spin ensemble are expressible solely in terms of symmetric processes 
and collective spin operators, the symmetric collective states of the ensemble are preserved. These 
many-body states, which are invariant under particle relabeling, can be efficiently simulated since 
they span a subspace whose dimension is linear in the number of spins. However, many open system 
dynamics break this symmetry, most notably when ensemble members undergo identical, but local, 
decoherence. In this paper, we extend the definition of symmetric collective states of an ensemble 
of spin-1/2 particles in order to efficiently describe these more general collective processes. The 
corresponding collective states span a subspace which grows quadratically with the number of spins. 
We also derive explicit formulae for expressing arbitrary identical, local decoherence in terms of 
these states. 
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I. INTRODUCTION 

The ability to model the open system dynamics of large 
spin ensembles is crucial to experiments that make use 
of many atoms, as is often the case in precision metrol- 
ogy GL S IS 0], quantum information science [!, 0, 0] 
and quantum optical simulations of condensed matter 
phenomena [1, 0, [l(| ■ Unfortunately, the mathematical 
description of large atomic spin systems is complicated 
by the fact that the dimension of the Hilbert space 
grows exponentially in the number of atoms TV. Realis- 
tic simulations of experiments quickly become intractable 
even for atom numbers smaller than N ~ 10. Current 
experiments, however, often work with atom numbers of 
more than N ~ 10 10 , meaning that direct simulation of 
these systems is well beyond feasible. Moreover, simu- 
lations over a range N ~ 1 — 10 are far from adequate 
to discern even the qualitative behavior that would be 
expected in the N 3> 1 limit. Fortunately, it is often 
the case that experiments involving large spin ensem- 
bles respect one or more dynamical symmetries that can 
be exploited to reduce the effective dimension of the en- 
semble's Hilbert space. One can then hope to achieve 
a sufficiently realistic model of experiments without an 
exponentially large description of the system. 

In particular, previous work has focused on the sym- 
metric collective states |"0s)i which are invariant un- 
der the permutation of particle labels: T\ij\ips) = IV's)- 
These states span the subspace M's C J^n, which grows 
linearly with the number of particles, dim(j£s) = Nj + 1. 
However, in order for Jf?s to be an invariant subspace, 
the dynamics of the system must be expressible solely 
in terms of symmetric processes, which are particle per- 
mutation invariant, and collective operators, which re- 
spect the irreducible representation structure of rota- 



tions on the spin ensemble. Fortunately, even within 
this restrictive class, a wide variety of phenomenon may 
be observed, including spin-squeezing 

turn and 

zero- 
temperature phase transitions [8j]. 

In practice, symmetric atomic dynamics are achieved 
by ensuring that there is identical coupling between all 
the atoms in the ensemble and the electromagnetic fields 
(optical, magnetic, microwave, etc.) used to both drive 
and observe the system (l3| . This approximation can 
be quite good for all of the coherent dynamics, because 
with sufficient laboratory effort, electromagnetic inten- 
sities can be made homogeneous, ensuring that interac- 
tions do not distinguish between different atoms in the 
ensemble. However, incoherent dynamics are often be- 
yond the experimenter's control. Although most types of 
decoherence are symmetric, they are not generally writ- 
ten using collective operators. Instead they are expressed 
as identical Lindblad operators for each spin 
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The fact that decoherence does not preserve J#s has 
been well appreciated and the standard practice in ex- 
periments that address the collective state of atomic en- 
sembles has been either: (i) to model such experiments 
only in a very short-time limit where decoherence can be 
approximately ignored; or (ii) to use decoherence mod- 
els that do respect the particle symmetry, but which are 
written using only collective operators. 

In this paper, we focus on ensembles of spin-1/2 parti- 
cles and generalize symmetric collective states to simply 
collective states, \ipc): which span the space J^c- This 
subspace, for which .^Cg C fflc C ^tv, is invariant un- 
der collective processes, which include symmetric deco- 
herence of the form of Eq. [1] We also give formulae for 
expressing arbitrary symmetric Lindblad operators in the 
collective state basis. Since dim(J4?c) = 0(N 2 ), this al- 
lows for efficient simulation of a broader class of collective 
spin dynamics and in particular, allows one to consider 
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the effects of decoherence on previous simulations of sym- 
metric collective spin states (see Ref. |l4|). We note that 
dynamics symmetries for spin- 1/2 particles have been 
studied in the context of decoherence- free quantum infor- 
mation processing (l5l . HH ] . Unlike our work, which uses 
symmetries to find a reduced description of a quantum 
system, these works seek to protect quantum information 
from decoherence by encoding within the degeneracies in- 
troduced by dynamical symmetries. 

The remainder of the paper is organized as follows. 
Section El reviews the representation theory of the rota- 
tion group, which plays an important role in defining the 
symmetries related to M's and 3%c- Section [ITT] intro- 
duces collective states and Section |IV] defines collective 
processes over these states. Section IVl gives an identity 
for expressing arbitrary symmetric Lindblads, e.g. Eq. 
[1] over the collective states. Section [VTl concludes. 



II. GENERAL STATES OF THE ENSEMBLE 

Consider an ensemble of N spin-1/2 particles, with the 
n th spin characterized by its angular momentum j' n ) = 
, ji n) , ji n) } ■ States of the spin ensemble are elements 
of the composite Hilbert space 

J4? N = ® JT (2) <g> • • • <g> Jt° {N) (2) 

with dim(^Ar) = 2 N . Pure states of the ensemble, \ip) G 
J%r, are written as 



w= E 



-mi .m2,...,mjv 



|mi,m 2 , ...,tojv) (3) 



with m n = ±~ and where 



|mi,m 2 , . . . ,m N ) = |i,mi}i(g)|i,m 2 }2®- • -® |^, m N ) N 

(4 ) 

satisfies 

j^\m 1 ,m 2 , . . . ,m N ) = hm n \mi,m 2 , ■ ■ .,m N ). (5) 

When studying the open-system dynamics of the spin en- 
semble, one must generally consider the density operator 



E 



m l i rn 2 > ■ -.jWiV 



x |mi,m 2 , . . . ,m N ){m' v m' 2 , . . .,m' N \ (6) 

States expanded as in Eqs. [5]and[6]are said to be written 
in the product basis. 



The basis kets | \ , m) for this particle therefore transform 
under the rotation R according to 



R\L m') = J2®i>,m( R )\l> m ) 



where the matrices *3 2 (R) have the elements 
@i, im = {l,rn'\R(a,0, 7 )\~,rn)- 



(8) 



(9) 



The rotation matrices ^2 (i?) form a 2— dimensional rep- 
resentation of the rotation group. 

For the ensemble of N spin-1/2 particles, each com- 
ponent of the ket \ip) = |mi, m 2 , . . . , tojv) transforms 
separately under a rotation so that an arbitrary state 
transforms as 



\^') = [^(R)f N \^). 



(10) 



The rotation matrices S$(R) — [@i(R)]® N provide a re- 
ducible representation for the rotation group but can be 
decomposed into irreducible representations (irreps) as 



■/max <4 
J=J min t=l 



(11) 



The quantum number i(J) = 1,2,..., d J N is used to dis- 
tinguish between the 



AH(2J+1) 



b N — JW 



(£ -J)!(f + J + 1)! 



, < J < 



N 
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degenerate irreps with total angular momentum J [1 71 ] - 
That is to say, d J N is the number of ways one can combine 
N spin-1/2 particles to obtain total angular momentum 
J. The matrix elements of a given irrep 2> J ' l (R) 

®m,m'(R) = (J,M,i\^(R)® N \J,M',i) (13) 

are written in terms of the total angular momentum 
eigenstates 



J 2 |J,M,i) = J(J + 1)| J, M, i) 
X\J,M,i) = M\J,M,i) 

with J z = J2n=i J» > J max = § and 



\ N odd 
N even 



(14) 
(15) 



(16) 



A. Representations of the Rotation Group 

For a single spin-1/2 particle, a spatial rotation 
through the Euler angles R = (a,/3,7) is described by 
the rotation operator 



It is important to note that degenerate irreps have iden- 
tical matrix elements, i.e. 

(J,M,i\9*{R) 9N \J,M',i) 

= {J,M,i'\^{R)® N \J,M',i l ) (17) 



R{a,(3,~f) = e 



-iajz -i/3j y -i-yjz 



(7) foralK,i'. 
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In this representation, pure states are written as 

-'max J d 'N 

M = E E J2c J>M!i \J,M,i) (18) 

J=,/ min M=-J i=l 

and mixed states as 

•/max -V 4,4' 

/5= E EE 

J,J'=J min M,M'=-J,J' i,i'=l 

pj,M,HJ>,M>,i>\J,M,i)(J',M , ,i'\ (19) 

States written in the form of Eqs. [TH] or [TO] are said to 
be written in the irrep basis. We stress that both the 
product and irrep bases can describe any arbitrary state 
in J1?n. 



III. COLLECTIVE STATES 

While the representations in Section |TT] allow us to ex- 
press any state of the ensemble of spin-1/2 particles, the 
irrep basis suggests a scenario in which we could restrict 
attention to a much smaller subspace of Jtfjsf. In particu- 
lar, the irrep structure of the rotation group, as expressed 
in Eq. [TTJ indicates that rotations on the ensemble do 
not mix irreps and that degenerate irreps transform iden- 
tically under a rotation. 

Following this line of reasoning, we introduce the col- 
lective states, \tpc)> which span the sub-Hilbert space 
J%c C J%n- Collective states have the property that 
degenerate irreps are identical; for pure states, cj t M,i = 
cj.m,i' for all i and i' . We note that the symmetric col- 
lective states mentioned in the introduction are the col- 
lective states with cj,m,i = unless J — and thus 
correspond to the largest J value irrep. We also note 
that 



dim Jifr 



J i(AT + 3)(iV+l), if ^ odd 
~ 1 ±(A + 2) 2 , if N even 



(20) 



Physically, the collective states reflect an inability to 
address different degenerate irreps of the same total J. 
This new symmetry allows us to effectively ignore the 
quantum number i and write 



Jmax J d N 

i^>= E E Ecj,m,*ij,m,z) 

J=J min M=-J i=l 

Jmax J 

= E E y d Ncj,M,N\J,M) ( 
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J=J min M=-J 

where we have defined effective basis kets 

dl 



with effective amplitude cj^m,n — cj,M,i fo r a ^ * (since 
the c.j y M,i are equal for collective states). 

At first glance, the factor of iJ d J N may seem bizarre. 
However, its presence enables us to apply standard spin- 
J operators to the effective kets without having to ex- 
plicitly refer to their constituent degenerate irrep kets 
\J,M,i). In other words, \ J,M,N) actually represents 
d J N degenerate kets, each with identical probability am- 
plitude coefficients. But since the matrix elements of 
degenerate irreps of a given J are identical, we need not 
evaluate them individually. 

As an example, consider a rotation operator R which 
necessarily respects the irrep structure of the rotation 
group. Calculating the expectation value of R by ex- 
panding the collective state \ific) m the irrep basis, we 
have 

= E E J2 c JM^J',M',AJ,M,i\R\J',M',i') (23) 

JJ' M,M' 

= E E J2 c J,M,i c JM'AJ,M,i\R\J,M',i) (24) 

J M,M> i 

= E E d N c *J,M,iW,i{J, M, 1\R\J, M', 1) (25) 
J M,M' 

where in going from Eq. [23] to [24] we set J = J' 
and i — i' since rotation group elements do not mix 
irreps. In reaching Eq. 1251 we have further used the 
collective state property that cj t M,i — cj,M',i'Vi, i' and 
the rotation irrep property that (J, M,i\R\J, M' ,i) = 
(J,M,i'\R\J,M',i')Vi,i'. We arbitrarily chose i = 1 as 
the representative element. 
Equivalently, if we define 

(J, M, N\R\J, M', N) = (J, M, l\R\J, M', 1) (26) 

and use the effective amplitude definition, we can 
evaluate the expectation using the effective basis kets 
\J,M,N) directly: 

<V>c|%c) 

= E E d J N c* J>M , N cj, M >,N(J,M,N\R\J',M',N) 

J. J' M.M' 

(27) 

= E E d J N cl MiN cj, M ,, N {J,M,N\R\J,M,N) . 



J M.M' 



(28) 



d J N factor allows us to calculate 



'd J N cj,M,N\J, M i N ) '= E C 'h M ,i\J, M, i) (22) 



We see then that the 

expectations directly in the collective notation. 

We can similarly define collective state density opera- 
tors, pc, which have the properties that (i) there are no 
coherences between different irrep blocks and (ii) degen- 
erate irrep blocks have identical density matrix elements. 
The second assumption again means we can effectively 
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drop the index i, since pj,M,i-,J,M',i ~ PJ,M,i';J,M',i' for 
any i and i' . This allows us to write 

Jmax J 

pc= E E 

Pj,m,n ; j,M',n\J, M, N)(J, M\ N\ (29) 

where the generalization of Eq. [22] for effective density 
matrix elements is given by 

d J N pj,M,N;j,M'-N\J, M, N)(J, M', N\ := 

X)pj,af, ( jj 1 jtf',i|J,M,t)<J > M / ,t| . (30) 

i=l 

Just as for the effective kets, the factor of d J N ensures 
expectations are correctly calculated using the standard 
spin- J operators. The density matrix has ^(N + 3) (AT + 
2)(N + 1) elements. 

We stress that Eq. [3DJ is different than naively taking 
the outer product of the effective kets defined in Eq. [22] 
That approach would involve outer products of kets be- 
tween different, although degenerate, irreps. Such terms 
are strictly forbidden by the first property of collective 
state density operators. Instead, one should consider 
both the effective kets and effective density operators as 
representing d J N identical copies of a spin- J particle. Eqs. 
|2"21 and [3DJ are then just rules for relating pure or mixed 
states of this effective spin- J particle to the true states 
in the irrep basis. 



IV. COLLECTIVE PROCESSES 

We are now interested in describing quantum pro- 
cesses, C, which preserve collective states, p' c = Cpc- 
Writing this explicitly, we must have 

E E ^p' JuMuN . JlMiiN \Ji,M 1 ,N)(J 1 ,M[,N\ 

Ji Mi,M( 

= E E d NPJ2M 2 ,N-J 2 Mi,NC\J2,M 2 ,N)(J 2 ,M^N\ 



J 2 M 2 ,Mo 



(31) 



If we define the action of C on collective density matrix 
elements as 



J,M,M',N 



C\J,M,N){J,M',N\ (32) 



we immediately see that this action must be expressible 
as 



/ 



J,M,M',N 



E E Xj^^uM^NHJuM^Nl 

Ji M 1 ,M[ 

(33) 

in order for the equality in Eq. [ST] to be met. Here 
^■"j^Mx M' IS an arbitrary function of its indices. Any 



process which preserves collective states by satisfying Eq. 
1331 is a collective process. 

Examples of collective processes are those involv- 
ing collective angular momentum operators { J x , J y , . . .} 
and more generally, arbitrary collective operators C = 
J2n=i ■ Since collective operators correspond to pre- 
cisely the rotations considered when defining the irrep 
structure of the rotation group, they can all be written 
as 

C = E E c J,m,m'\Ji M, N)(J, M', N\ , (34) 

J M,M' 

which cannot couple effective matrix elements with dif- 
ferent J. 

However, the collective operators define a more restric- 
tive class than an arbitrary collective process, which can 
couple different J blocks, so long as it does not create 
coherences between them. In fact, if all operators are col- 
lective, then the symmetric collective states (\ips)) span 
an invariant subspace of the map. This holds even when 
considering Lindblad operators that are written in terms 
of collective operators, such as collective spontaneous em- 
mission 



J-pJ+ - \ j+J-P - \pJ+J- 



(35) 



which describes a process in which all spins of the ensem- 
ble emit together. 

In the following section, we demonstrate that a process 
of the form 



JV 



■J,M,M' ,N 



^2a^\J,M,N)(J,M',N\(a^ , (36) 



which cannot be written solely in terms of collective op- 
erators, is nonetheless a collective process. This result 
was inspired by considering symmetric spontaneous em- 
mission 



Eii^if-^ 



i 



(n)3(n) 



J- P~ 2-7+ J 



(37) 

which corresponds to the spins emitting independently 
via identical processes. Accordingly, we call processes 
which are invariant under exchanging particle labels sym- 
metric processes. 

In order to appreciate how £r sym is related to Eq. [3DJ 
expand the single spin operator a in the spherical basis 



a = aj 



i + E a iii = ai ^ + a + J + + a -3- + azJz ( 38 ) 



1) A /oo" 
,00/' j~ — \10, 



with the convention h — 1, j+ 
and j z = ( J -i ) ■ We can write an arbitrary symmetric 
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Lindblad for a, of which Cr aym is a specific example, as 



n=l 
N 

= E 



a(n)p(a(n))t 



1 i . 1.1 

2 A NP - 2^ N 



(39) 



with the collective operator An given by 

N 

a n = Y,fr (n W n) 



= (-|o- 



(40) 



+ (a*_ai — a*_a z + et}a+ + a*a+)j+ 
+ (aj<2_ + a^a/ + a^a^ — a*a_) j_ 

+ (^l a -| 2 - ^l a +! 2 + a 7 a * + a *z a i)Jz ■ 

In this form, it is clear that only the first term of the 
symmetric Lindbladian is not written using collective op- 
erators. In fact, if we again expand a (n) in the spherical 
basis, we observe that the only terms which involve non- 
collective operators are those which do not involve the 
identity operator, 



N r 



E 

n=l 



a(«)^(a< n ))t 



= \ai\ 2 Np + ^2(a q a*jJ q p + a/a*/5J^) 

N r 

+ E E^^ 1 ) 1 



n=l 



q.r 



(41) 



N = 1 



1x2 



1 x 



lxl 



3x1 




lxf 



2x § 




1 x 



2x0 



FIG. 1: Degeneracy structure from adding spin-1/2 particles, 
labeled as dir x J 



with j q = 2nLii« a collective spin operator. In the 
following section, we demonstrate that these processes 
preserve the collective states and we give an identity for 
evaluating such processes for arbitrary spherical basis el- 
ements. 



V. IDENTITY 

Identity 1. Given a collective density matrix element 
for N spin-1/2 particles, \J,M,N)(J,M',N\, we have 



jJ,M,M',N 
N 

= J2ji n) \J,M,N)(J,M',N\0W 

n=l 
1 

~2J 



1 



Y J+1 
N 



2J+1 



d J N J+l 



A J q > M \J,M q ,N)(J,M^N\A 



(42) 
(43) 

J,M' 



+-^j B i' M \ J ~ L M «> N )( J - !. <> N \ B r M ' 



y J+1 



<2(J 



—D J q M \J+ 1, M g , N)(J + 1, M' r , N\D J r 



J,M' 



(44) 



where q,r £ {+, -, z}, M + = M + 1. M_ = M - 1 and 
M, = M, 



*N = E d N = 



AT! 



j'=j 



:i-j)!(f 



and 



J, A/ 



7,M 



V(J- M)(J + Af + 1) 
V(J + Af)(J-Af + 1) 
M 



Af)(J-M-l) 
-V(J + M)(J + Af-l) 



73/' M = + M)( J -M) 



and lastly 



£>+' M = - v / (J + M+l)(J + M + 2) 
J ' M = y/(J - M +1)(J - M + 2) 



73 
73 



7,M 



= V(J + M + l)(J-Af + l) 



(45) 



(46a) 
(46b) 
(46c) 

(47a) 
(47b) 
(47c) 

(48a) 
(48b) 
(48c) 



Note that a N and d^ are zero if J is negative or J = 
tV/2, ensuring that only valid density matrix elements 
are involved. 

In the following subsections, we prove Identity [1] in- 
ductively. The motivation for the inductive proof comes 
from the simple recursive structure of adding spin-1/2 
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particles. As seen in Fig. [D the d J N irreps which corre- 
spond to a total spin J particle composed of TV spin- 1/2 
particles can be split into two groups, depending on how 
angular momentum was added to reach them. By ex- 
pressing the N particle states in terms of bipartite states 
of a single spin-1/2 particle and a spin-(JV — 1) particle, 
we can then evaluate the dynamics independently on ei- 
ther half by assuming Identity [T] holds. Returning the 
resulting state to the N particle basis should then con- 
firm the Identity. By inspection, the base case of N = 1 
holds, as the A^ M terms reduce to the single spin-1/2 
matrix elements. We now proceed to the inductive case. 

A. Recursive state structure 

In order to apply the inductive hypothesis, we need to 
express an N particle state in terms of N — 1 particle 
states. This recursive structure is best seen by examin- 
ing Fig. [TJ which illustrates the branching structure for 
adding spin-1/2 particles. For example, the three- fold de- 
generate N = 4 spin-1 irreps arise from two different spin 
additions — adding a single spin-1/2 particle to the non- 
degenerate J = |,iV = 3 irrep and adding to the 2-fold 



degenerate J = |, N = 3 irreps. Since we are always 
adding a spin-1/2 particle, the tree is at most binary. 
This allows us to recursively decompose the degenerate 
irreps for a given J in terms of adding a single spin-1/2 
particle to the two related N — 1 degenerate irreps. 

Recall that for the collective states, wc define effective 
density matrix elements which group degenerate irreps 
(Eq. |3"0|) . In order to make the relationship between 
states of different N more obvious, we rewrite the state 
\J,M,i) as \J,M,N,i). The N and i indices indicate 
the state is from i-th degenerate total spin-J irrep that 
comes from adding N spin-1/2 particles. So that we can 
leverage the binary branching structure seen in Fig. [TJ we 
need to then relate the N particle irrep states to the N—l 
particle irrep states. Accordingly, we define J, M; \ , J± 
\,N — where the last four entries indicate that 

the overall N spin state can be viewed as combining a 
single spin-1/2 particle with a spin J ± i particle. The 
spin J ± i particle is from the i\-st such irrep for N — l 
spin-1/2 particles. With these definitions, we can now 
relate the N particle states to the N—l particle states 
by explicitly tensoring out a single spin-1/2 particle: 



\J,M,N)(J,M',N\ 
1 4 

^^|J,A/,iV,z)(J,M',7V,z| 



-v i=1 

-''N — l 



v il= l 



N ia =l 

J+i 
N-l 



di 



N 



E 



J,M (jz> mi 



(49) 
(50) 



1 I I I I 

= rfT E \J,M;- 7 J+-,N-l 7 i 1 )(J,M';-,J+-,N-l,i 1 \ 



1 II I I 

-j- \J,M;-,J - -,N - l,i 2 )(J,M';-,J - ~,N - l,i 2 \ 



(51) 



E 



(i mi|<J + i M'-m' 1 ,JV-l| J ^(7*' mi 



J+i,M'-mi 



'I ,m' 2 \(J- -,M'-m' ,N -l\ J > M 'cl^\ f , , 

'2 2IN 9 J—^,M'—m' 2 



(52) 



with Clebsch-Gordan coefficients JM Cf^™\ 
(J, M;ji,j2\ji,mi;j2,rri2) and the rtii,m! i sums over 
single spin projection values ±i. In reaching Eq. 1521 we 
made use of the definition of the effective density matrix 
element for N — 1 spins given in Eq. With this 

recursive state definition, we can now start the inductive 
step of the proof. 

B. Applying inductive hypothesis 

In order to prove the Identity, we must be able to ap- 
ply the inductive hypothesis to Eq. Ignoring the 



Clesbsch-Gordan coefficients for the moment, consider 
an arbitrary term from Eq. 1521 The dynamics distribute 
as 
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A' 



(n) 



|i, m i )(^,m' i \®\J±^M-m i ,N-l)(J±^, M-m^iV-lj 



» 



=/ (g) |j± i,M - mi,N- 1}(J ± i,M- mj, JV - 1| + |~, m^XI ® f^-muM'-mlN-i (53) 



r 



By extension, all terms in Eq. 



split the dy- / 



J±&,M-mi,M'-m'i,N-l 



This means evaluating the / co- 



namics in this manner, which allows us to apply efficients according to Eq. 1441 after which we rewrite the 
the inductive hypothesis to evaluate f2' m i' m i^ and bipartite states in the N spin basis. 



We have the /i'" 11 '™!' 1 terms 
1 



d N-l J+l r 

— V V V A^' miJl < M »C^ miq J,M(]i' mi \L M-~ T+- N-l 

N j 1= l J 1 = J mi L 



x £ 2 [ m; ; i, J + i, N - 1, h\ J ' M 'Cf^ M , 

J[—J 

and the fi' m2 ^ m 2> 1 terms 



(54) 



j',,M'rih m ir |-'"' 
°J+i Af'-m' A 



w i 2 =l J 2 = J-1 m 2 
J 



«2 



2^ {J 2 ,M r ,-,J- -,N-l,l 2 \ ^j^M>- m ^ C J-±,M>-m> 2 Ar 
J' 2 —J—l m' 2 

The /J+i.M-mi,M'-mi,iV-i terms are 



(55) 



4(2J+1) 



j+ 3 

_^ a N— i 2,7 -)- 2 



■/+§ 

"iV-l 



E EE 

il=l Ji=J mi 

.7+1 , 
x V- V- IV *,r>. 



J[—J m[ 



r h m 'i a j - 

^J+h,ML-m\ Ar 



(56) 



7+i 



,J-, 



ii = l J\=J— 1 mi L 
J 



J+i.Af-mij M |,mi J.M-^l.mi T » . 1 7 1 at 1 ,• \ 



X 



— r 1 ™/ L 



J+i M'-l 



J{— J— 1 m^ 



(57) 



JV-1 



44 + j2(j+f) jl=1 Jl= 



%_1 J + 2 

E E E 



2) ii=l Ji=J+l mi 
J+2 

x E E 

J{— J+l 



^J+|,M-mij liM9 ^i,mi 



J+|,M'-m' 



J+i,Af'-t 



(58) 
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and lastly, the / 



J-k,M-m-2,M' -m'„,N-l 



terms are 



E E E 



J 

x E E 

J^ — J—l mi, 



A 



/' M CjT M _ m2 |J 2 ,M,;-,J--,7V-l,i 2 > 



2' 



1 ^' r' o' o' ' ■ a ' J — 4* ,M' —mL J— ± M'— m„ 



J-l.M' -m 



(59) 



J- 3 
l JV-l 



CjV-1 J_i 

/•' /''" 4 ->, 7 : , EE 

a N a N-l /, \ J 2/ »a=l Ja=J— 2 m 2 
J-l 

x E E 

J2 — J— 2 m~ 



9 J-|,M g -m2 J-i,M-m 2 



I J 2 ,M,;-, J --,JV-l,i 2 > 



; \J,M' rth m 2 



4,M^-,J-~,N-l,i 2 y>™C^? M/ . 



c 

J-f ,M'-m 



J— 4,M'— m 



(60) 



.7+1 

K JV-1 



<%-l J+l 



J_I_I z rf z a 

J+l 

x EE 

J' — J mi, 



E EE 



J-i ,M-m 2 J 2 t M„ fik'™ 2 



/~i 2 ' 



Cy^ M . m2 \J2,M q ,- J+ - N~l^ 2 ) 



\p M'- 1 71 1 Af 1 j \J,M'n3' m S J' 2 ,M' rr h m 2 n J-h,M-m' 2 

(61) 



r 



Evaluate sums 



We are now tasked with showing that Eqs. I54ll6ll sum 
to f J ' M ' M > N as written in Eq. 23] Before doing so, 
we observe that the Ji , rrii and J[ , m! i sums factor in all 
the equations above. Moreover, if one replaces primed 
quantities with unprimed ones, the Clebsch-Gordan and 
A, B, D coefficients of the kets in a given Ji, rrii sum are 
identical to those of the bras in the related J'i,m'i sum. 
Therefore, we focus on simplifying the unprimed sums 
and then apply those results to the primed sums in order 



to simplify Eqs. [5HI6T1 In Appendix [3] we explicitly cal- 
culate two representative sums from these equations. The 
calculations involve manipulating products of Clebsch- 
Gordan and A, B, D coefficients. Although tedious, the 
interested and pertinacious reader should have no trouble 
evaluating them for all relevant sums, finding in partic- 
ular that the J ± 2 terms vanish. We forego detailing 
all those manipulations here and simply use the results 
in both the primed and unprimed terms of the equations 
above, which then simplify 



J 



Eq. EH to 



1 11 11 

rf j (2J + 2)2 E D i M \ J + !> 2' J + 2 ' N ~ liH){J + hK; r J+ r N ~ h illD '' M 

~A J q M \J, M q ; J + 1 N - 1, h)(J + 1, M' r ; i J+ i JV - l,ii\D? M ' 
-D J q M \J + 1, M q ; l,J+l,N- Mi><J, M' r ; I J + i,N - 1, ii|A^' M ' 
+Af 7 | J, M g5 i J + i TV - 1, m; ; i J+ i JV - 1, n|^ M ' , 



(62) 
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Eq. [55] to 



"3 

"N-l 



I I I I I 

— J2 B i' M \ J - l ^J-pN-i, - 1, m; ; -, J — -; JV — 1, h\B J r > M ' 



+^ M | J, M q ; ~, J - i; JV - l.ixXJ - 1, m; ; I, J - i; JV - 1, h|S/' m ' (63) 
+B q ' M \J - 1, M 9 ; ~, J - ~; TV - 1, h)(J, M' r ; ~, J - ~; JV - 1, ?1 |^ M ' 
+A q ' M | J , Af g ; ~, J - 1; AT - l,i 1 )(J, J<; i, J — i; AT — 1, h|^ M ' , 



Eq. [56] to 



J + 3 



2(J+|) 



^<' M | J, M g ; i J + 1 AT - 1, i x )(J + 1,m; ; ~, J + \, JV - 1, ii|D r J ' M ' 



(2J + 2) 2 9 1 ' "'2' 2' ' ±/x ' r '2' 2 

2(J + ^) 11 11 

' J2J + W D < ' |J+ M<?; 2 ' J + 2 ' ^ ~ lj Zl)(J ' M ^ 2' 3 + 2' ^ " h lM "' M 
(J+l) 2 



——A J q ' M \J, M q ; 1,J+~,N-1, h){J, M' r - 1, J+l, N-l, i x \A J r ' M ' 



(64) 



Eq. [57] to 



J+i 

a N-l 



J— — 

d J N d N _\2J(J f ^ » 1= i 



11 11 

r £ (J + 1)B^ M | J - 1, M q ; -, J — -; JV — 1, h){J - l,M' r ; -,</—-; JV - 1, i \ B J r < M 



A J q M \J, M q ; \, J - \-N - 1, h)(J - 1, m; ; i J - i; JV - l,h\B^ M ' ( 65 ) 



2' 2 ^—n 2 :» 2 

+Bg' M |J - 1, M g ; ~, J - i; JV - 1, ii)(J, J<; ~, J - i; JV - 1, n] A r J 

7^Y^' M I 4 M q , i, j - i- at - 1, ^(J, m; ; i J - 1; JV - l,ii|A^ r 



Eq. [58] to (since J + 2 terms vanish) 



rv J+1 1 13 13 

<2( J +l)<+i ^ M |^ + l,M g ;-,J+-,JV-l,i 1 }(J+l,M; ; - ) J+-,JV-l ) i 1 | J D^ , 

Eq. [59] to 



(66) 



^Tj^yrry [i + £ b>"\j - i,m,\,j -In - 1,^-1,^^-^-1,^ 



-2(J - i)A^ M | J, M 9 ; 1 J - 1; JV - 1, i x >(J - 1, M' r \ i J - ±; JV - 1, 

-2(J - i)J3^ M | J - 1, M g ; 1 J - i; JV - 1, ^(J, J<; 1, J - I; JV - 1, i x \A J r ' M ' 

+4( J - i) 2 ^ M | J, M g ; I, J - 1; JV - LiiXJ, M;; 1 J - \; JV - 1, h\A^ M ' , 



(67) 
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Eq. [60] to (since J — 2 terms vanish) 

i 3 
d 5 

d * N 2Jd .,-i E BJ> M \J-l,M q] -,J--,N-lM){J-l,M>--,J--,N-lMB. 



N a N i 1= l 



and Eq. [61] to 

J+h 



,J+i 



7 

=1 


J + 


1 


+ 


1 






J + 


1 




1 




+ 







,J.M\ 



1 1 



1 

r 
i 



yJ,M> 

\J,M' 



A J q M \J, M q - ~ J + ~ jv - 1, n}( J, m; ; ~ J + ~ AT - 1, i^Ai'M' 



J(J + 1) 



(68) 



(69) 



D. Recover / J 



1 



1 



N+ 1 



< 8J 2 + 20J+12 



1 



We now combine the equations from the previous sub- 
section to recover the Identity in Eq. [44] Given that 
density operators in the collective state representation 
lack coherences between different J irreps, we expect 
| J ± 1)(J\ and \J)(J ± 1| terms to vanish. Since both 
the | J){J ± 1| and | J ± 1) ( J| terms have the same coeffi- 
cients, we need only explicitly deal with one of the two. 

Starting with \ J+1)(J\ coefficients from Eqs. [62][64]and which gives overall 
[691 we find 



d J N (2J + 2) 2 \ (2J + 3)(2J+1) 
1 2J + 7V + 4 



J+l 
N 



<2(j+i) 4+ 1 



J{N + 2J + 2) 
2J+ 1 



(74) 



1 



1 



(2J + 3) 



<V (2J + 2) 2 (2J+l)(2J + 2) : 



u JV-l 2 



*7V-1 



(J+l)(2J+l)^ + j 



1 



<(2J + 2) 2 

: 



1 (iV + 1) | N + 2J + 2\ 



2J+1 



2J+1 



Similarly, for |J — 1)(J| coefficients in Eqs. 
1571 we have 



1 



d J N ±.P 



i Q ^ 2J -[1- 



2 J 



Turning to J + 1 terms from Eqs 
coefficients sum to 



C ] J 2 J + i 

u jv-i 2 



J (70) 

(71) 
and 



= (72) 
andRM the 



1 



1 



d-'\(2J + 2) 2 (2J+ l)(2J + 2) : 



J_|_3 

a N— i 2«7 + 2-, 



2 

(73) 



J 



J+i 

*JV-1 



(J+1)(2J+ 1)^+4 



/+ 1 

N 



d J N 2(J + 1) 4+ 1 



1 11 

^ £ D^|J+l,M,;-,J+- J JV-l,i 1 ) 

(75) 



x (j + i,m; ; ~ j + 1 at -i,^!^ 



The J terms from Eqs. 



1 and l6"9l have coefficients 



+ 



(2J + 3) 2 



d J N \{2J + 2) 2 (2J+l)(2J + 2) ; 



■[1 



C^N— 1 2<7 ~h 2 -j 
,J+i J + 3 J 
"iV-1 2 



+ 



a 



N-l 



J(J+1)(2J+1) ^+5 

N — l 



d J N {2J + 2f 



(7V + l)(2J + 3) 7V + 2J + 2 



2J+ 1 



J(2J+ 1) 



d J N 2J 



4 +1 



2J+ 1 



< J+l 



(76) 
(77) 
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which gives overall 



1 

2J 



a J N +1 2 J + 1 
d J N J+l 



jj- £ A J q < M \J,M q ;\, J+l,N-l,h) 



i 2 =l 



x (J,M; ; i J+i JV-l,ii|^ M ' . (78) 

EU and [57] have 



Similarly, the J terms from Eqs. 
coefficients 



1 



J+k 



1 



2 J 



2 ? J + o 



"<2J 



2 J 



7V-2J 
2J+1 V + l 



1 



a 



J+i 

N 



2J+ 1 



d J N J+l 



which gives 



1 

2J 



/ +1 

N 



2J+ 1 



< J+l 



2J-1 



(79) 



From the definition of \J, M, N){J, M' , N\ given in 
Eq. [3D1 we see that Eqs. [75J and [50] correspond to the 
| J, M, N)(J, M', N\ terms in Eq. [44] A similar combina- 
tion of Eqs. [551 and 1521 corresponds to the J— 1 term and 
the combination of Eqs. [66j and [75] corresponds to the 
J term. We have thus shown inductively that Identity Q] 
holds. ■ 



VI. CONCLUSION 

We have presented an exact formula for efficiently ex- 
pressing symmetric processes of an ensemble of spin-1/2 
particles. The efficiency is achieved by generalizing the 
notion of collective spin states to be any such state which 
does not distinguish degenerate irreps. For a collection 
of N spin-1/2 particles, the effective Hilbert space di- 
mension grows as N 2 , a drastic reduction from the full 
Hilbert space scaling of 2 N . The collective representation 
is used in Identity [U which gives a closed- form expres- 
sion for evaluating non-collective terms from symmetric 
Lindblad operators. Unfortunately, due to the compli- 
cated structure of adding spin- J > i particles [l7|, our 
results do not appear to generalize. Nonetheless, we be- 
lieve that this approach will become a useful tool in ana- 
lyzing collective spin phenomenon and in particular, ac- 
curately considering the role of decoherence in collective 
spin experiments [14| . 



1 1 1 

-j- A J q > M \J,M q ;-,J--,N-l,h) 



d J 



x (J,M^, J-i N-l,h\Ai> M ' 



(80) 



And finally, the J — 1 sums from Eqs. [UJ [55] and 1571 

have coefficients 



<4P 



\2J{J+l) 



+ 



1 



2 ( J - 2 

1 



X N-1 



2 J 



/-I J 

1 



N -2J 



2J-1 2J+1 



J+l 



1 



1 



N 



dpJd^ 1 



2J- 1' 
(81) 



which gives 

^ £ B J q > M \J-l,M q ;\,J-\, N-l,h) 

(82) 



d J N 2Jd J N 



i 1 = l 



x (J-l,M; ; i J-±,N-l,h\B J r > M ' 



APPENDIX A: EXPLICIT SIMPLIFICATION OF 
TYPICAL SUMS 



In IV CI we simplify the sums in Eqs. 15415T1 but do 
not go through the detailed algebra. The work involves 
manipulating products of Clebsch-Gordan and A,B,D 
coefficients. In this appendix, we explicitly calculate two 
representative sums from this set and invite the reader 
to calculate the remainder in a similar fashion. 

First, consider the sums over J\ and mi in Eq. 1541 
which is representative of sums in Eqs. [54] and [55] For 
Ji = J + 1 



A; 



<C 



aI'-^+^cV, , 



2(J + 1) 



J,M (~<2'~2 

-y/(J + M + 2)(J + M +1) q= + 

q = - 

q = z 



y/ (J- M + 2)(J- M + 1) 

y/(J-M+l)(J + M+l) 



z _t_ D J,M 

2J + 2 q 



(Al) 
(A2) 
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and for ,J\ = J 



Similarly, for J\ — J, we have 



j+\,M-h 



T I , 1 1 

J,MQ2'2 



J+hM-h 



A\ 



■J.M, 



q (-l2> 2 q J.M 



1 



2(J+1) 

1 



y /(J-M)(J + M+l) q 
y/(J + M)(J-M+l) q 
M q 



A>> M 

2J + 2 9 



i i i 
where A\' 2 = A< 



(A3) 



(A4) 



= 0, Ai 



= A 2 J 2 =1 and 



Similarly, consider the sums over J\ and m\ in Eq. [57j 
which is representative of Eqs. For J\ = J — 1, we 

have 



B 



2 ' 



J, M £J 2 ' 2 



J,M^f2' _ 2 



'(J-Af g )(J-M + l) B , 7+ i,M-i 



4J(J + 1) 




J + M 



J + M+ lBg + i,M+ ^ 
J - M V J — M + 1 b ^+5.m-5 



\j-M q ){J- M + 1) D j+i,M-i 2(J + 1) 



4J(J + 1) 



Bo 



J - M + 1 



'J+l 



y/( J-Af)(J-M-l) g = + 
V(J + M)(J + M - 1) g=- 
v/(J + Af)(J-M) g = z 



+ 1 TJ.J.M 



(A5) 



B, 



J+i,M-| JM, 



11 

J,MQ2'2 



iC 2 ' 2 



' (J + Af g )(J-M + l) g ,/ + i,M-i 



4J(J+ 1) 



J-M 9 J + M + 1B, 



j+M q y j - m + i B 



J+hM-k 



{J + M q ){J - M + 1) b j+i,m-1 



4J(J+ 1) 
l 




M±3j 



+ 



V /(J-M)(J + M + 1) g = + 
^(J + M)(J- M + 1) q=- 
M q = z 



J(J+l) Aq 



J,M 



(A6) 
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